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Abstract 

In analyzing the nature of thermal radiance experienced by an accelerated observer 
(Unruh effect), an eternal black hole (Hawking effect) and in certain types of cosmolog- 
ical expansion, one of us proposed a unifying viewpoint that these can be understood 
as arising from the vacuum fluctuations of the quantum field being subjected to an 
exponential scale transformation in these systems. This viewpoint, together with our 
recently developed stochastic theory of particle-field interaction understood as quantum 
open systems described by the influence functional formalism, can be used effectively to 
address situations where the spacetime possesses an event horizon only asymptotically, 
or none at all. Examples studied here include detectors moving at uniform acceleration 
only asymptotically or for a finite time, a moving mirror, and a two-dimensional col- 
lapsing mass. We show that in such systems radiance indeed is observed, albeit not in 
a precise Planckian spectrum. The deviation therefrom is determined by a parameter 
which measures the departure from uniform acceleration or from exact exponential ex- 
pansion. These results are expected to be useful for the investigation of non-equilibrium 
black hole thermodynamics and the linear response regime of backreaction problems 
in semiclassical gravity. 
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1 Introduction 



Particle production |j] with a thermal spectrum from black holes [g, |3], f|], moving mirrors 
0, accelerated detectors |J, observers in de Sitter Universe [0] and certain cosmological 
spacetimes || has been a subject of continual discussion since the mid-seventies because 
of its extraordinary nature and its basic theoretical value. The mainstream approach to 
these problems relied on thermodyamic arguments |9|, |10[, finite temperature field theory 
techniques [|TT| , 
spacetime) 



13 1 , or geometric constructions (event horizon as a global property of 



FJ], or pairwise combinations thereof. The status of work on quantum field the- 

The eighties saw attempts and 



ory in curved spacetimes up to 1980 can be found in 
preparations for the backreaction problem [TJJ (for cosmological backreaction problems, see 
flTf), i.e., the calculation of the energy momentum tensor (see |18| and earlier references), the 
effect of particle creation on a black hole (in a box, to ensure quasi-equilibrium with its radi- 
and the dynamical origin of black hole entropy [19]. These inquiries are mainly 



ation) H 

confined to equilibrium thermodynamics or finite-temperature field theory conditions.^] To 
treat problems of a dynamical nature such as the backreaction of Hawking radiation on 
black hole collapse, one needs a new conceptual framework and a more powerful formal- 
ism for tackling non-equilibrium conditions and high energy (trans-Planckian) processes. A 
new viewpoint which stresses the local, kinematic nature of these processes rather than the 

1, which regards the 



FT 



traditional global geometric properties has been proposed 
Hawking-Unruh thermal radiance observed in one vacuum as resulting from exponential red- 
shifting of quantum noise of another. This view puts the nature of thermal radiance in the 
two classes of spacetimes on the same footing [ 3^] , and enpowers one to tackle situations 
which do not possess an event horizon at all, as the examples in this and a companian paper 
will show. 

Such a formalism of statistical field theory has been developed by one of us and co- 
workers in recent years |T5], [T7], | 



|39|j . This approach aims to provide the quantum 
statistical underpinnings of field theory in curved spacetime, and strives at a microscopic 
and elemental description of the structure and dynamics of matter and spacetime. The start- 
ing point is the quantum and thermal fluctuations for fields, the focus is on the evolution 
of the reduced density matrix of an open system (or the equivalent distribution or Wigner 
functions); the quantities of interest are the noise and dissipation kernels contained in the 
influence functional |4(|, and the equation of motion takes the form of a master, Langevin, 
Fokker-Planck, or stochastic Schrodinger equation describing the evolution of the quantum 
statistical state of the system, including, in addition to the quantum field effects like ra- 
diative corrections and renormalization, also statistical dynamical effects like decoherence, 
correlation and dissipation. Since it contains the causal (Schwinger-Keldysh) effective action 
f4Tfl it is a generalization of the traditional scheme of thermal field theory |13| and the 'in-out' 
(Schwinger-DeWitt) effective action [|42]| , and is particularly suited for treating fluctuations 
and dissipation in backreaction problems in semiclassical gravity PU . 



The foundation of this approach has been constructed recently based on the open system 



1 Among other notable alternatives, we'd like to mention Sciama's dissipative system approach |p0f , Un- 
ruh's work on sonic black holes |2l| (see also Jacobso n p2[ ), Zurek and Thome's degree of freedom counts 
p3[, Sorkin's geometric or 'entanglement' entropy pj 



theory approach 



(see also |26|]), and Bekenstein-Page's information 
See also the views expressed earlier by Stephens, t'Hooft and Whiting p9j. 
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concepts and the quantum Brownian model E3, |3 



The method has since been applied 
to particle creation and backreaction processes in cosmological spacetimes [[15|, [|(| [47|. For 
particle creation in spacetimes with event horizons, such as for an accelerated observer and 
black holes, this method derives the Hawking and Unruh effect [48] , |44j from the viewpoint 
of exponential amplification of quantum noise |33[] . It can also describe the linear response 
regime of backreaction viewed as a fluctuation-dissipation relation. |47], [49| . 

This paper is a continuation of our earlier work [48], [|9| to present two main points: 

1) A unified approach to treat thermal particle creation from both spacetimes with and with- 
out event horizons [34] based on the interpretation that the thermal radiance can be viewed 
as resulting from quantum noise of the field being amplified by an exponential scale trans- 
formation in these systems (in specific vacuum states) [Q. In contradistinction to viewing 
these as global, geometric effects, this viewpoint emphasizes the kinematic effect of scaling 
on the vacuum in altering the relative weight of quantum versus thermal fluctuations. 

2) An approximation scheme to show that near-thermal radiation is emitted from systems un- 
dergoing near-uniform acceleration or in slightly perturbed spacetimes. We wish to demon- 
strate the relative ease in constructing perturbation theory using the statistical field theory 
methods. Let us elaborate on these two points somewhat. 

It may appear that this approximation can be equally implemented by taking the conven- 
tional viewpoints (notably the geometric viewpoint), and the perturbative calculation can be 
performed by other existing methods (notably the thermal field theory method). But as we 
will show here, it is not as easy as it appears. Conceptually, the geometric viewpoint assumes 
that a sufficient condition for the appearance of Hawking radiation is the existence of an 
event horizon, which is considered as a global property of the spacetime or the system. (Note 
for the case of an extreme Reissner-Nordstrom spacetime, this is not the there exists 

an event horizon but no radiation. Q) When the spacetime deviates from the eternal black 
hole, or that the trajectory deviates from the uniformly accelerated one, physical reasoning 
tells us that the Hawking or Unruh radiation should still exist, albeit with a non-thermal 
spectrum. But the event horizon, if exists, of the deformed spacetime may not be so easily 
described in geometric terms. And for time-dependent perturbations of lesser symmetry, or 
for situations where uniform acceleration occurs only for a finite interval of time, it is not 
easy to deduce the form of Hawking radiation in terms of purely global geometric quantities 
(see, however, Wald ]50| and Teitelboim et al [|5l|). The concept of an approximate event 
horizon, which exists for a finite period of time or only asymptotically, is difficult to define 
and even if it is possible (by apparent horizons, e.g., f52|), rather unwieldy to implement in 



the calculation of particle creation and backreaction effects. 

Technically one may think calculations via the thermal field theory are equally possi- 
ble. Indeed this has been tried before by one of us and others. One way is to assume a 
quasi-periodic condition on the Green function, making it near-thermal [O. But this is not 



a good solution, as the deviation from eternal black hole or uniform acceleration disables 
the Euclidean section in the spacetimes (Kruskal or Rindler), and the imaginary time fi- 
nite temperature theory is not well defined any more. Besides, to deal with the statistical 



2 This is the one case which arose in the discussion between Hu and Unruh (private communication, 1988) 
who shared the somewhat unconventional view that the exponential redshifting is a more basic mechanism 
than event horizons responsible for thermal radiance. 
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dynamics of the system, one should use an in-in boundary condition and work with causal 
Green functions. The lesson we learned in treating the backreaction problems of particle 
creation in cosmological spacetimes [j^, |45| is that one can no longer rely on methods which 
are restricted to equilibrium conditions (like the imaginary-time or thermo-field dynamics 
methods), but should use nonequilibrium methods such as the Schwinger-Keldysh (closed 
time path) effective action [^TJ for the treatment of backreactions. Its close equivalent, the 
influence functional method |^0| is most appropriate for investigating the statistical me- 
chanical aspects of matter and geometry, like the entropy of quantum fields and spacetimes, 
information flow, coherence loss, etc. |53| . 



In this and a companian paper [|54j] we shall use these methods to analyze particle cre- 
ation in perturbed situations whose background spacetime possess an event horizon, such as 
an asymptotically uniformly-accelerated detector, or one accelerated in a finite time interval 
(Sec. 3), the moving mirror (Sec. 4) and the asymptotically Schwarzschild spacetime (Sec. 
5). In the follow-up paper we shall study near-thermal particle creation in an exponen- 
tially expanding universe, a slow-roll inflationary universe, and a universe in asymptotically- 
exponential expansion. What ties the problem of thermal radiation in cosmological as well as 
black hole spacetimes together is the exponential scale-transformation viewpoint expressed 
earlier BoL M, M M M. 



The stochastic theory approach is capable of implementing this 
view. One can describe all these systems with a single parameter measuring the deviation 
from uniformity or stationarity, and show that the same parameter also appears in the near- 
thermal behavior of particle creation in all these systems. This result will be used later in 
our exploration of the linear-response regime of the backreaction problem in semiclassical 
gravity. 



1.1 General Formalism 

Consider a particle detector linearly coupled to a quantum field. The dynamics of the 
internal coordinate Q of the detector in a wide class of spacetimes is derived in ||44fl , and can 
be described by Langevin equations of the form : 

fir d r)J ft 

3Q-Jt^- 2 L/ (Md ^ m (L1) 

where £(£) is a stochastic force with correlator (£(£)£(£')) = hu{t,t'). The trajectory x M (s) 
of the detector, parametrized by a suitable parameter s, will be denoted simply by x(s) for 
convenience. 

For the special case of linear coupling between a field <ft an d the detector of the form 
Lint = gQ4>(x(s)), the kernels /x and z/, called the dissipation and noise kernels respectively, 
are given by 

2 2 

l*M = jG(x(s),x(s'))^- l j<[^x(s))J(x(s'))]> (1.2) 

u(s,s') = e lG^\x(s),x(s'))^j<{^x(s))J(x(s'))}> (1.3) 

where G and G^ are the Schwinger and the Hadamard functions of the free field operator 
evaluated for two points on the detector trajectory, < > denotes expectation value with 
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respect to a vacuum state at some arbitrarily chosen initial time U, and [ , ] and { , } denote 
the commutator and anticommutator respectively. This result may be obtained either by 
integrating out the field degrees of freedom as in the Feynman- Vernon influence functional 
approach [40] or via manipulations of the coupled detector-field Heisenberg equations of 
motion in the canonical operator approach. 

It will often be convenient to express the kernels /i and v as the real and imaginary parts 
of a complex kernel £ = v + i/j,, called the influence kernel. For linear couplings, it follows 
from the above expressions that ( is given by the Wightman function G + : 

C(s, s') = e 2 G + (x(s),x(s')) = e 2 < > • (1-4) 

The influence kernel thus admits the mode function representation 

as,s') = e 2 Y,M^s))ul(x(s')), (1.5) 

k 

the w^s being the fundamental modes satisfying the field equations and defining the partic- 
ular Fock space whose vacuum state is the one chosen above. This method of evaluating the 
kernels /i and v is only applicable for linear coupling cases. 

An alternative approach [44], consists of decomposing the field Lagrangian into para- 
metric oscillator Lagrangians at the very outset, thus converting a quantum field-theoretic 
problem to a quantum mechanical one. Denoting the parametric oscillator degrees of free- 
dom by q k (and their masses and frequencies by m k and uj k respectively), the detector-field 
interaction mentioned earlier is generally given by L int = J2k c k(s)Qq k , where the coupling 
"constants" c k now become time-dependent, and contain information about the detector 
trajectory. In this approach, the influence kernel is given in terms of the oscillator mode 
functions X k , as 

C(s, s') = f°° dk I(k, s, s')X k ( S )X* k (s') (1.6) 
Jo 

where the X k s satisfy the parametric oscillator equations 

X k + uftt)X k = (1.7) 

satisfying the initial conditions X k (ti) = 1 and X k (ti) = —iu)(ti). The spectral density 
function I(k, s, s') is defined as 

77 >\ ^ \ c k (s)c k (s') 

I(w,s,s) =}^5{uj-uj k )- (1.8) 

V 2m k (t i )uj k (t i ) 

One may decompose the influence kernel into its real and imaginary parts, thus obtaining 
the dissipation and noise kernels: 

1 r°° 

l*M = * dkI(k,s, S ')[X* k ( S )X k ( S ')-X* k ( S ')X k (s)} (1.9) 

2 Jo 

1 r°° 

u(s,s') = - dkI(k,s,s')[X* k (s)X k (s')+X* k (s')X k (s)]. (1.10) 

2 Jo 

By expressing the field as a collection of parametric oscillators, it can be explicitly verified 
that the two approaches mentioned above lead to the same result for the influence kernel £. 
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For the purpose of calculating it in a specific case, we will find it more convenient to use the 
second approach. 

To study the thermal properties of the radiation measured by a detector, the influence 
kernel is compared to that of a thermal bath of static oscillators each in a coherent state 



0: 



poo 

(= / dkI eS {k,E) [C k {Z)coskA-ismkA] (1.11) 
Jo 

where 

£ = (t + t')/2, A = t-t' 

and the function C k = coth -^f- We will show in the specific cases discussed below that 
the unknown function C k indeed has a coth form, and can then deduce the temperature of 
the radiation seen by the detector. Here I e s(k, E) is the effective spectral density, also to be 
determined by formal manipulations of (1.6). We can always write ( in this way since v is 
even in A while \x is odd. By equating the real and imaginary parts of the two forms of ( 
and fourier inverting, we obtain 



oc 



I cS C k = -I dA cos£;A u(E, A) (1.12) 

7T J-oo 



'off 



/oo 
dA sin£;A //(£, A). (1.13) 
-oo 



I f°° 
71 



We will now consider various examples where ( is evaluated and shown to have, to zeroth 
order, a thermal form. Higher-order corrections to ( give a near-thermal spectrum. In princi- 
ple, the real and imaginary parts of the influence kernel may be substituted in the Langevin 
equation (1.1) to yield stochastic near-thermal fluctuations of the detector coordinate Q. 
This procedure will be demonstrated in the example of a finite time uniformly accelerated 
detector (Sec. 3 below). In this way, the methodology presented above describes a stochastic 
approach to the problem of detector response, as opposed to the usual perturbation theory 
approach (where the perturbation parameter is e 2 ) involving the calculation of detector tran- 
sition probabilities. It should be emphasized that equation (1.1) is exact for linear coupling 
and does not involve a perturbation expansion in e 2 (for linear systems, such an expansion 
is, in fact, unnecessary because they are exactly solvable). 
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2 Asymptotically Uniformly- Accelerated Observer 



We first consider the case of a non-uniformly accelerated monopole detector in 1 + 1 dimen- 
sions. For a general detector trajectory (x(r),t(r)) parametrized by the proper time r, it 
has been shown [J{| that the function ((t, t') is 



roc rife 

£( r T >) = v + ifJL = |_ / ^ e - ifcWT) - t(r ' )) cosA;(a;(r)-x(r / )). (2.1) 

Here e is the coupling constant of the detector to a massless scalar field (initially in its 
ground state). The initial state of the detector is unspecified at the moment and would 
appear as a boundary condition on the equation of motion of the detector. Here, however, 
we are primarily interested in the noise and dissipation kernels themselves, as properties of 
the field, and not in the state of the detector. 

First, we note that the function ( can be separated into advanced and retarded parts, in 
terms of the advanced and retarded null coordinates v(r) = t(r)+x(r) and u(r) = t(r)—x(r) 
respectively: 

C(r,r') = f-J^e- ik ^- v ^ (2.2) 



Ati Jo k 
C(r,r') = C(r,T')+CXr,r'). 

In the case when the detector is uniformly accelerated with acceleration a, its trajectory 
is given by: 

v(t) = -e aT ; u{r) = -e~ aT . (2.3) 
a a 

Substitution of the above trajectory into equations (2.2) yields a thermal, isotropic detector 
response at the Unruh temperature o/(27r) |44, f49| . 



2.1 Perturbation increasing with time 

The above analysis is now applied to the case of near-uniform acceleration by introducing a 
dimensionless h parameter which measures the departure from exact uniform acceleration: 

h = - 2.4 

cr 

where the overdot indicates derivative with respect to the proper time. The trajectory of 
the detector is now chosen to be: 

v(t) = -^-ef a{T)dT ; u (t) = -^—e- $ a{T)dT . (2.5) 
a(r) a{T) 

One can expand a(r) in a Taylor series about the acceleration at r = 0: 



<r) = a + E -l a o ] (2-6) 



n=i n - 
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where a n ^ denotes the ra-th derivative of a at r = 0. We make the assumption of ignoring 
second and higher derivatives of a. This implies 



o(r) = ao + h ral 



(2.7) 



where h = a'o/a 2 ,. 

Hereafter, we shall also make the further assumption of evaluating quantities to first order 
in ho- In this approximation, h — ho to first order in ho- Then there is no distinction between 
h and ho (h is essentially constant), and we can safely drop the subscript and work with 
h alone. It should be noted that expanding quantities to first order in h actually involves 
expansion of quantities to first order in hrao, and hence, for arbitrary trajectories, the final 
results are to be considered valid over time scales r such that r <C (/ia ) _1 . Alternatively, 
equation (2.7) can be taken to define a family of trajectories for which this analysis applies. 

Using the linearized form of a(r), one can now obtain the trajectory explicitly, to first 
order in h. The result is: 



v T 



u r 



V 0T ( 1 + hTCLoi— ~ 1 



(1 hra (^f- ■ 1) 



(2.8) 



One also finds, to first order in h, 



-ik(v(r)— t'(r')) 



a V 2 I 



1 _ ikhe ^ ( (^! + aoS 2 _ 2S) sinh( ^ ; 



+ A(a S - 1) cosh( 



a A, 



-ik(u(r)— u(t')) 



e a o 2 ' 



A(a £ + 1) cosh( 



2 

1 + ikhe 
a A 



(2.9) 



- aoS ((^ + a S 2 + 2S)sinh(^) 



(2.10) 



where A = r — r', S = |(r + r') 
Using the identities H8j 



e a o 



^ e - a os sinh( ^) 4 f°° 



,2k 



dvK2iv( — e a ° s )[cosh(7ri/) cos(z/aoA) — 2sinh(7ri/) sin(i/aoA)] 
it Jo ao 

(2.11) 



and 



T(iu) 



7T i^/l \ ■ o 7T 

. r ; I r(- + | 2 = — - — 

usmtntP 2 cosli7rz/ 



(2.12) 



one finally obtains, after some simplification 

e 2 f°° dk r . ,irk 



( a (r,r / ) = — / ^[coth(— )cos(fcA)(l + / i r 1 ) -isin(JfeA)] (2.13) 
47r jo k ao 

p2 /.qo A h. Trie 

?( T y) = ?L — [coth(— )cos(fcA)(l + hr 2 ) -zsin(fcA)] (2.14) 

47T JO k an 
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with 



2 7r£\ 



T l = -£;tan(£;A)tanh 2 ( — 

a 



/c tan(ZcA) tanh 



2/ ttA;, 
a ' 



, a A 2 
' 4 
a A 2 



,a A 



a A, 



+ a S 2 - 2E) sinh(^— ) + A(a S - 1) cosh(^— ) 



+ a S 2 + 2E) sinh( 



a A, 



A(a E + l)cosh(^) 



(2.15) 



The advanced and retarded parts of Re(Q being unequal, the noise is anisotropic. Adding 
expressions (2.12) and (2.13), we have 



POO fjU TTfC 

C(r,r') = — / ^[coth(— )cos(fcA)(l + fcr) -isin(fcA)] (2.16) 
2n Jo k an 



rk 
a 



where 



ri + r 2 



,irk 



a A 



a A, 



/cStan(ZcA) tanh ( — )(2 sinh — a A cosh — — ) 



(2.17) 



2 v ' v a 2 2 

The noise experienced by the detector is thus identical to the noise experienced in a heat 
bath, with a small correction, T. The accelerated detector therefore has a near-thermal 
response at the Unruh temperature a /(2n) with an order h correction which increases with 
time. 



2.2 Perturbation exponentially decreasing with time 

We will now consider a trajectory for the accelerated detector which exponentially approaches 
the uniformly accelerated trajectory at late times. This trajectory, in null coordinates, has 
the form 

v (t) = 1 e aoT (l + ae~ lr ); u{t) = -a^e'^il + ae~ lr ). (2.18) 
In this case, the magnitude of the proper acceleration is, to first order in a, 



?y2 

a(r) = a {l + mT 7T (l + — )} + 0{a 2 ). (2.19) 



The influence kernel is obtained in a manner similar to the treatment of the previous sub- 
section. Here, we get 

e 2 r°° dk r . ,nk. 



( a (r,T , ) = — / ^[coth(— )cos(fcA)(l + ar 1 )-2sin(A;A)] (2.20) 
ATT Jo k do 

p 1 ^ roo pfV- ttJc 

?( T y) = — / ^[coth(— )cos(A;A)(l + ar 2 )-2sin(A;A)] (2.21) 
4% Jo k an 



with 



T, = -2ka^e-^ sinh (a ° 7) A tan(fcA) tanh 2 (— ) 

2 do 

T 2 = -2ka V 7 ^ sinh (q ° + 7) A tan(fcA) tanh 2 (— ) (2.22) 

2 a 

9 



The noise is again seen to be anisotropic. Adding ( a and ( r , we have 

2 r°° dk r , , 7T k , 

a 

where 



C(r,r') = - ^[coth( — )cos(fcA)(l + ar)-^sin(A;A)] (2.23) 
2ir Jo k an 



T = £i±ll = _2A;a - 1 e- 7s sinh ^ cosh ^ tan(ifcA) tanh 2 (^). (2.24) 

In this case, the correction to the thermal spectrum is exponentially suppressed at late times. 
This feature will distinguish the behavior of quantum fields in the vicinity of a moving mirror 
and a collapsing mass, as shown in later sections. 
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3 Finite-Time Uniformly- Accelerated Detector 



In this section, we consider a detector trajectory which is a uniformly accelerated one for a 
finite interval of time (— to, to). Before and after this interval, the trajectory is taken to be 
inertial, at uniform velocity. To ensure continuity of the proper time along this trajectory, 
the velocity of the detector is assumed to vary continuously at the junctions ±£o- 
With these constraints, the trajectory is chosen to be 

x(t) = XQ 1 (a~ 2 — tot), t < —to 

= (a~ 2 + t 2 )i t > -t , t < to 

= x^ia- 2 + t t), t>t . (3.1) 

The trajectory is symmetric under the interchange t — > —t. a is the magnitude of the 
proper acceleration during the uniformly accelerated interval (—to, to) of Minkowski time 
and Xq is the position of the detector at time t . x and t are related by Xq — t^ = a~ 2 . 
Before the uniformly accelerated interval, the detector has a uniform velocity — t /x (we 
have chosen units such that c = 1; if one keeps factors of c, the velocity is —c 2 to/x ) and 
after this interval, its velocity is to/x . This trajectory thus describes an observer traveling 
at constant velocity, then turning around and traveling with the same speed in the opposite 
direction. The "turn-around" interval corresponds to the interval of uniform acceleration. 

We may also define null coordinates u = t — x and v — t + x. In terms of these, the time 
at which the trajectory crosses the future horizon u — of the uniformly accelerated interval 
is tn = — (a 2 iio) -1 - 

If we choose to parametrize the trajectory by the proper time r, it can be expressed as 
(with the zero of proper time chosen at t — 0) 

u(t) = 9(-T - r)v {a(r + r ) - 1} - a" 1 ^ + t)#(t - r) e - aT 

+6(t - to)u {1 + a(r - r)} (3.2) 

v(r) = -9(-T - r)u {a(r + r ) + 1} + a'^To + r)9(r - r)e ar 

+6(t - roMl - a(r - r)} (3.3) 

where ±r is the proper time of the trajectory when it exits (enters) the uniformly accelerated 
phase. It satisfies the relations 

v = t + xo = a~ l e ar ° 

uo = t - x = -a~V aT °. (3.4) 

Another convenient definition is the horizon - crossing proper time ±t# = ±(a~ x + r ). 

The function ((r, r') can be found in a standard way. If both points lie on the inertial 
sector of the trajectory, it has the usual zero-temperature form in two-dimensional Minkowski 
space. If both points lie on the uniformly accelerated sector, it has a finite temperature form 
exhibiting the Unruh temperature. It is therefore straightforward to obtain the following: 
If r, t' > t or r, r' < — r , 
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If -Tq < T,T' < Tp , 

£(r, t') = f°° ^{coth(— ) cos k(r' - t) + i sin k(r' - r)}. (3.6) 
27T jo k a 

Also, if r < — t , t' > To, 

POO rjV- 

(( T , t') = — ^ cos fc((r' + r) tanh(ar ))e lfc(r - r+2(a tanh(aro) - ro)) . (3.7) 

Z7T JO 

Of interest is this function evaluated for one point on the inertial sector and the other on the 
uniformly accelerated sector. We will show that this function has a thermal form if the point 
on the inertial sector is sufficiently close to (£0,^0) and departs smoothly from the thermal 
form away from it. It is also found that the horizons of the uniformly accelerated sector 
(which are not horizons for the entire trajectory) are the points where the near-thermal 
expansion breaks down. 

Consider, for example, the case when — rp < t' < t and r < — r . Then the function ( is 
expressed as 

£( T) r ') = f_ J°° ^^-ikia^e-^' +vo(a(T+T )-l)) + ^Ha^e^' +u a (a(r+r a )+l)) | 

Introducing the Fourier transforms 

e* e = - — / due luJT F( )(-)a e 2 0) k > 

2na J-00 a a 



^T(-)(-)^ e ff, A;>0 (3.9) 
J-00 a a 



27ra j- 

we obtain, after some simplification, 

C(r,r') = f- {cos k (t' + t + - \n(l - a(r + t ))) coth(— ) 

Att Joky V a / a 

+ 2 sin (V + r + ^ln(l - a(r + r ))^ 

1 / nk \ 

+ cos fc(r' + r In ] a(r + th) |) coth( — )0(r H + r) + 0(-t h - r) 

a \ a / 

+ % sin fc(r' + r — - In I a(r + t h ) \)9(t~h + t) \ . (3.10) 
a J 

If we further restrict our attention to the case r > —th, i.e. both points lie inside the Rindler 
wedge, the above expression simplifies to the following : 

e 2 r°° dk ( k 



«^'> - i?rf cos u in(i - a2(T+T ° )2) ) x 

{71k { 1 
coth( — ) cos k r' + T H ln(l - 2 
a V 2a 



T + Tq 
T + T H > 

r + r H 



+ *sin k (t' + t + J- 111(1-2^1))}. (3.11) 
V 2a t + Tu J J 
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It is clear from this expression that an exact thermal spectrum is recovered in the limit of 
r — > — t , as expected. Suppose we now define r + t = e as the time difference between the 
proper time r and the proper time of entry into the accelerated phase, — tq. Then ae will 
be the appropriate dimensionless parameter characterizing a near-thermal expansion. Note 
that e < 0. 

From the above expression for (, we find that there is no correction to the thermal form 
of C(T)T') to first order in e. This can be understood from the fact that the coordinate 
difference between the point r — — t — e and a corresponding point on a globally uniformly 
accelerated trajectory with the same proper time is of order e 2 . Indeed, we may define 
Rindler coordinates {ip,v) 011 the right Rindler wedge by v = ^~ 1 e 5 ™ and u = — ■0~ 1 e~"^"'. 
Then the Rindler coordinates for the point t = — t — e on the trajectory we consider are 
found to be ip = a + 0(e 2 ) and r\ = — t — e + 0(e 2 ), which are exactly the coordinates, 
to order e, of a corresponding point with the same proper time on a globally uniformly 
accelerated trajectory with acceleration a. It is thus no surprise that the spectrum is exactly 
thermal up to order e. 

Furthermore, it can be shown in a straightforward way from the above expression that the 
spectrum is also thermal up to 0(e 2 ), although the above-mentioned coordinate difference 
does have terms of order e 2 . Then the first correction to the thermal spectrum is of order e 3 
and has the form 

C(t,t') = ^ / ^{coth(-)cosA;(T'-T + ^f )-*sin£;(T'-T+^)} + 0(e 4 ). (3.12) 
271 Jo k a A 6 

The validity of such a near-thermal expansion is characterized by the requirement that 
| ae | is small. This translates to — 1 < a(r + To) or equivalently, r > —th- The expansion 
thus breaks down for r < —th, for which case the two-point function may be called strictly 
non-thermal. This is the case when one of the points lies outside the right Rindler wedge 
while the other point is still inside it. The two-point function in such a situation will contain 
non-trivial correlations across the Rindler horizon, as was pointed out before f49|. 

The response of the detector is governed by the Langevin equation (1.1). This equation 
may be formally integrated to yield 

(Q(t)Q(t')} = ^ f ds f rfsV(s,s , )e- 7(T - s) e- 7(T '- s ' ) sinn(T-s)sinfi(T / -s / ) (3.13) 

i ' J — OO J —CO 

where f2 = (Qq — 7 2 )^, fio is the natural frequency of the internal detector coordinate and 
7 = e 2 /4 is the dissipation constant arising out of the detector's coupling to the field. The 
double integral in the above equation may be computed by splitting each integral into a part 
which lies completely in the uniformly accelerated sector and parts which lie in the inertial 
sectors. For example, suppose we wish to compute the above correlation function for the 
case —to < t, t' < To, i.e. both points lie in the uniformly accelerated sector. Then each 
integral can be split into two parts (JL^ = Jl^ + jT T0 ) and the resulting double integral 
therefore has four terms: 

(Q(r)Q(r')} = F 1 + F 2 + F 3 + F 4 . (3.14) 
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Writing v = Re((), we obtain, after straightforward manipulations, 

f 1 = ^> r° ds r ds > r ^-^-^-^ x 

7TJ6 J — oo J — oo JO rv 

sin fi(r — s) sin f2(r' — s') 

= ^_ e -7(r+r' + 2r ) ^ [( 2 _ p + Q 2 )2 + ^2]-! x 

nir Jo k 

{(tt 2 + 7 2 + k 2 ) cos Q(r - t') + (fi 2 - 7 2 - k 2 ) cos fi(r + r' + 2r ) 
+2 7 fisinfi(r + r +2r )}, (3.15) 



and 



F 4 ^itef d s r'^r^e^'-)coth(^)e-(-)e-^'-') 
7ri 2 J-To J-T Jo k a 



x 



sinfi(r — s) sinf2(r' — s') 

= i^_ e -7(r + r' + 2. () ) ^ dk r (Q 2 + 2 _ £.2)2 + ^2]-! x 

nil 2 Jo k a 

{( 7 2 + A; 2 + fi 2 ) cos Q(r - r') + (fi 2 - 7 2 - A; 2 ) cos Q(r + r' + 2r ) 

+27^(sinfi(r + r' + 2r ) - sinfi(r + r ) - sinfi(r' + r )) 

-fi 2 (cosfi(r + r ) + cosfi(r' + r ) - 1)}, (3.16) 

where Re stands for the real part. 

The functions F2 and F3, in which one of the integration variables runs over the inertial 
sector and the other over the uniformly accelerated sector, are difficult to evaluate. We shall 
simply express them here in the following form : 



— TO 

dse 13 sin ft(r-s) x 

-00 



POO A lr 

/ r ik(au s+uo(l+aT )) A , k .n + e -ik(av s-v (l-aT )) A J k . /g-^X 

Jo k 

F 3 = ^L e -^ T+T '^Re f T ° dse^ s sin n(r' - s) x 

7TJ I J — oo 

f°° ^\ e -ik(*u s+u (l+aT )) A r k . T \ + ^(at^-t^l-aTi,)) r )] /g^g) 
Jo ft 

where the functions Ai and A 2 are 

Ai(ife;s) = f ds , e ifco_le ' M 'e 78 'sinn(s-s / ) 

A 2 {k;s) = f ds'e- lka ~ le ~ as ' ' e /s ' sinfi(s - s'). (3.19) 

J— TO 

Similarly, if one wishes to compute the detector correlation function for two points in the 
late inertial sector (r, t' > tq) , then one has nine terms similar in form to the ones displayed 
above. 
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4 Moving mirror in Minkowski space 



In this section, we treat the motion of a mirror following a trajectory z(t) in Minkowski 
space. A massless scalar field is coupled to the mirror via a reflection boundary condition. 
It obeys the Klein-Gordon equation 



subject to the boundary condition 



^-^ = (4.1) 

dt 2 dx 2 y ' 



4>(t,z{t)) = 0. (4.2) 



For a general mirror path this equation is difficult to solve; however we can exploit the 
invariance of the wave equation under a conformal transformation to change to simpler 
coordinates. We follow the treatment of H. To this end, we introduce a transformation 
between the null coordinates u, v and u, v defined as 

u = t — X, V = t + X, 

u = f{u), v = v (4.3) 

The function / is chosen such that the mirror trajectory is mapped to z = 0. To do this, we 
relate the two sets of coordinates as follows: 

* = l[v + f(u)} 

x = \[v-f{u)] (4.4) 
On the mirror path, setting z = means that the trajectory can be expressed as 

\[t-f(t)] = z(l[t + f(t)]) (4.5) 

which allows / to be implicitly determined. In the new coordinates the wave equation is 
unchanged, however it now has a time independent boundary condition, meaning the mirror 
is static, while the detector moves along some more complicated path. Thus the wave 
equation with boundary condition can easily be solved to give 

(f)(t,x)= (2nky 1/2 smkx e~ lkt dk (4.6) 
Jo 

where the mode functions are orthonormal in the Klein-Gordon inner product. In these 
barred coordinates, £ is proportional to the two point function in the presence of a static 
reflecting boundary at z — 0. 

Also, in these coordinates, the time dependent modes of the field are just exponentials. 
That is, they can be described by oscillators with unit mass and frequency k. So Xk(t) is 
a solution to the oscillator equation (1.4), and by satisfying the initial conditions Xk(0) = 
l,X' h (0) = —ik we obtain 

X k (t) = e- m (4.7) 



15 



We now consider a detector placed in the vicinity of the mirror. The spectral density 
function / is determined by the path of the detector and its coupling to the field. Denoting the 
detector position by r(t) and the field modes by qk(t) and assuming the monopole interaction 



J 'mt 



eQ4>(t, x) 5(r — x) dx 
eQ4>(t, r) 

J eQlk(t) sinkrdk, (4i 



we have 



I(k,t,t') = ( ^^-5{k — k n )e 2 sin k n r{t) sin k n r{t' 



1 II 

e 2 

= — sin kr(t) sin kr(t ). (4-9) 
Defining u = t — f(t) and v = t + r, we can now express the function ( as 

£ = e2 I °° ^ Lik{u'-u) _ e ik{u'-v) _ e ik(v'-u) _|_ e ife(v'-«)l (4-10) 



7T JO 

Since only the outgoing modes have reflected off the mirror, only the outgoing part of the 
correlations ( will give appropriate thermal behavior. Thus, from now on, we focus on the 
correlation 

u=— r d ±^'-^. (4.H) 

Sti Jo k 

It remains to evaluate the above function. To do this, we specify the function / by considering 
a specific mirror trajectory. A convenient choice of the mirror path is the following: 

z (t) = -t- Ae' 2Kt + B (4.12) 

for A, B, k positive. This path possesses a future horizon in the sense that there is a last 
ingoing ray which the mirror will reflect; all later rays never catch up with the mirror and so 
are not reflected. It is this aspect which enables the moving mirror to emulate a black hole. 
Eq. C |4=.5f) can now be solved to give 

f(t) = -t-Un^. (4.13) 

In the late time limit (t ~ B), we consider the following ansatz for f" 1 

f-\x) ~ B- Ae~ K{B+x) + a (4.14) 

where a is taken to be small in the sense that terms of order a 2 are ignored. In this 
approximation, one finds 

a = - K A 2 e~ 2K{ - B+x) (4.15) 
and the transformation from barred to unbarred coordinates becomes 



u 



B - Ae~ K{B+u) - KA 2 e- 2K{B+u) (4.16) 
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plus terms of higher powers in e - K ( B + u ). 

We now need an explicit form for the detector trajectory u(t) since this is what appears 
in the function (. Choosing it to be inertial, we have r(t) = r* + wt, which gives u(t) = 
t(l —w) — r*. In terms of the proper time of the detector, this becomes u(r) = t^J (|^) — t*- 

Defining the sum and difference E = |(r + r') and A = t — r', and 2; = ^ we 
obtain 

^ _ u = _2Ae- K(B+r * +&) sinh( — ) - 2nA 2 e- 2 < B+r * + ^ sinh(«i/A). (4.17) 
This is substituted in and, after some simplification we obtain the near-thermal form 

p2 /»oo AT™ Trie 

Cuu(r,r') = — / — [coth(— ) cos(fcA)(l + hV) -ism(kA)} (4.18) 

o7T JO fc f€Z/ 

with 

T = -2ke^ B+r ' +J]v) tanh 2 (— ) tanfcA. (4.19) 

Thus a thermal detector response, at the temperature Doppler-shifted by a factor z/ 
depending on the speed of the detector, is observed, with a correction that exponentially 
decays to zero at late times. 
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5 Collapsing Mass in Two Dimensions 



In this section we study radiance from a collapsing mass, analogous to the moving mirror 
model. We essentially follow the method of , but using stochastic analysis, and generaliz- 
ing it to include higher order terms in the various Taylor expansions involved, thus exhibiting 
the near-thermal properties of detector response. 

We will exploit the conformal flatness of two dimensional spacetime in the subsequent 
analysis. Outside the body the metric is expressed as 

ds 2 Q = C(r)dudv (5.1) 

where u, v are the null coordinates 



a 



t-r* + R* (5.2) 



v = t + r* + R* 
and r* is the Regge- Wheeler coordinate: 

r dr' 



C(r') 



(5.3) 



with Rq being a constant. The metric outside the body is thus assumed to be static in order 
to mimic the four dimensional spherically symmetric case (for which Birkhoff's theorem 
holds). The point at which the conformal factor C = represents the horizon, and the 
asymptotic flatness condition is imposed by C — > 1 as r — > oo. 

On the other hand, the metric inside the ball is for now assumed to be a completely 
general conformally flat metric: 

dsl = A(U, V)dUdV (5.4) 

with 

U = T-r + R 

V = T-r-R (5.5) 

and Ro and Rq are related in the same way as r and r*. The surface of the collapsing ball 
will be taken to follow the worldline r = R(r), such that, for r < 0, R(r) = Rq. Thus, at 
the onset of collapse, r = t = 0, U = V = u = v = 0on the surface of the ball. 

We will let the two sets of coordinates be related by the transformation equations 

U = a(u) 

v = f3(V). (5.6) 

The functions a and f3 are not independent of each other because one coordinate transfor- 
mation has already been specified by equation ( |5.3| ). 

Without as yet determining the precise form of a and (3, we will consider a massless 
scalar field propagating in this spacetime subject to a reflection condition 0(r = 0, r) = 0. 
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Such a field propagates in a similar fashion to the field in the vicinity of a moving mirror. 
To make this explicit, we introduce a new set of barred coordinates 



u = (3[a(u) - 2R ] 

v — v. (5.7) 

In terms of these, we also define the coordinates r — \{y — u), and t — |(U + u). 

These new coordinates have the properties : a) r = =>- r = 0, and b) the field equations 
have incoming mode solutions of the form e tkv . Thus the left-moving parts of the correlation 
functions of the 'in' vacuum defined in terms of barred coordinates are identical to those of 
the vacuum defined with respect to unbarred coordinates. 

Keeping these properties in mind, we may expand the field in terms of standard modes 
obeying the reflection boundary condition (by conformal invariance of the massless scalar 
field equation) as 

<!>(r,t) = Jj'£ 1 lk(t)s™kr, (5.8) 

v L k>0 

just as in the moving mirror case. 

We now consider a detector placed outside the collapsing ball at fixed r (or r*), namely 
r = r (or r* = rj). The interaction between detector and field is described by the interaction 
Lagrangian 

Lint = -eQ(j)(s,r), (5.9) 

where 

r = ^(v-/3(a(u)-2R )) (5.10) 

= ±(t + r* -R* -P(a(t-r* + R* )-2Ro)) 

s = ^(t + r*-R* + (3(a(t-r* + R* )-2R )) 

and Q is the internal detector coordinate. 

The influence kernel (, due to a reflection condition at r — 0, has the same form as the 
moving mirror case, in barred coordinates. Its outgoing part is therefore given by 

?^-> (5.n) 

where 

u = s - r = P(a(t - r* -R* )- 2R ) (5.12) 

and u' is the same function of t' . 

We will now determine the functions a and (3 and show that, to zeroth order in an 
appropriate parameter, u is an exponential function of t, and thus ( uu has a thermal form. 
The correction to the exponential form, obtained by including higher order terms, will lead 
to a near-thermal spectrum. 
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To determine a and (3 we match the interior and exterior metrics at the collapsing surface 
r = R(t). Then we have 

= -C^^[l + (l + ^(l-fl 2 ))*]-' (5.13) 
(fV) = % = C-'^H-d + ^d-fl 2 ))*] (5.14) 

where R = 4^- 

Now we expand these quantities about the horizon. We recall the definition of the horizon 
radius Rh as C(Rh) = 0. We may further define r h as R{r h ) = R h . Then we obtain the 
following Taylor expansions: 

R(t) = R h + u(r h - r) + (3(r h - rf + • • • (5.15) 

where v = —Rfa), [3 = \Rijh), and 

f)C 1 r) 2 C 

= 2Kv(T h -T) + (2Kp + 1 v 2 )(r h -T) 2 + ••■ 

where k = the surface gravity, and 7 = Since the ball is collapsing, v > 0. 

Substituting the above expansions in the expression for a'(tt), we obtain, to order (r^— r) 2 , 

^ = a(R -R h + Th -U)+ b(R -R h + r h - Uf (5.17) 
du 

where 

a = (p + 1)k (5.18) 

6 = 1{(3 + u)P+(1 + u)^-Uk(1-u 2 )(1 + u)}. (5.19) 
i> 2 

Note that, for a slowly collapsing ball, !/<l, and hence a reduces to the surface gravity k. 
Also, to order {jh — r), 

/77J 

— = c + d(r h + R h -R - V) (5.20) 



where 

A(l + v) 



2v 

A Ak 



(5.21) 



d = ^(^- T (l-^)(l + 4 (5-22) 

We consider a regime in which {r^ — r)d <C c so that we may ignore the second term in 
( |5.20| ). Then we can integrate this equation to give 

v(V) = /3(V) = ci + cV (5.23) 
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where c\ is an integration constant. 

Similarly to lowest order in \ (which turns out to be the appropriate dimensionless 
parameter describing deviations from exact exponential scaling or exact thermal behavior), 
we integrate equation ( [5.171) to give 



U{u) = a(u) = R -R h + T h + a~ l e-< u - c "\l + -^e" a(n - C2) ), (5.24) 



c 2 being another integration constant. 

We are now in a position to obtain explicitly the transformation between barred and 
unbarred coordinates, to lowest order in X. Thus we have 

u = (3[a{u) - 2Rq] 

= Mi + M 2 e~ a{u - C2) {I + -^e" a( "- C2) ) (5.25) 



a? 



where 



Mt = Cl - c(Ro + R h - r h ) (5.26) 

M 2 = -. (5.27) 
d 

At the position Tq of the detector, u = t — rg. Therefore, defining A = u' — u and S = 
|(m' + u) + Tq, we may perform the above transformation to obtain 

u'-u = -2M 2 e aC2 {e- a(s - r S ) sinh ^- + A e -2a(s-r 5 -f ) sinhaA }, ( 5 2 8) 



Invoking the identities ( |2.11| ) and ( p,12|) , the function ( uu can now be simplified to yield the 
near-thermal form 

Cuu = — ^{coth( — )cosA;A(l + T) -isinkA} (5.29) 
m Jo k a 

where 

9hk irk 

r = -^e a(c2_E_r o ) tanh 2 (— ) tan A;A sinh aA. (5.30) 
a 6 a 

The function V vanishes at late times (X — > oo). Thus the exact thermal spectrum is 
recovered at the Hawking temperature redshifted by the velocity of the surface of the ball, 
on a time scale defined by the surface gravity a. 
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6 Discussion 



We now summarize our findings and discuss their implications. There are four main points 
made or illustrated here: 



1) This paper gives a stochastic theoretical derivation of particle creation, in the class 
of spacetimes which possess an event horizon in some limit. This approach generalizes the 
established methods of quantum field theory and thermal field theory (in curved spacetimes) 
to statistical and stochastic field theory. The exact thermal radiance cases arising from an 
exact exponential scale transformation such as is found in a uniformly-accelerated detector, 
the Schwarzschild black hole and the de Sitter universe, have been treated in the stochastic 
theoretical method before ||44| . Here we give the treatment of the moving mirror and the 
collapsing mass as further examples. (Thermal radiation in certain classes of cosmological 
spacetimes inflationary universe will be studied in a following paper.) 



2) We have shown that in all the examples considered in this class of spacetimes, i.e., 
accelerated observers, moving mirrors and collapsing masses (black holes), those which yield 
a thermal spectrum of created particles all involve an exponential scale transformation. Ther- 
mal radiance observed in one vacuum arises from the exponential scaling of the quantum 
fluctuations (noise) in another vacuum. This view espoused by one of us |St] 
illustrated in the examples treated here. 
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3) The main point of this paper is to show how one can calculate particle creation in 
the near-exponential cases, yielding near-thermal spectra. These cases are not so easy to 
formulate conceptually using the traditional methods: the geometric picture in terms of the 
properties of the event horizons as global geometric entities works well for equilibrium ther- 
modynamics (actually thermostatics) conditions, so does thermal field theory which assumes 
a priori a finite temperature condition (e.g., periodic boundary condition on the imaginary 
time). They cannot be generalized to non-equilibrium dynamical conditions so easily. In the 
stochastic theory approach we used, the starting point is the vacuum fluctuations of quan- 
tum fields subjected to kinematical or dynamical excitations. There is no explicit use of the 
global geometric properties of spacetimes: the event horizons are generated kinematically 
by exponential scaling. (Thus, for example, this method can describe the situations where a 
detector is accelerated only for a short duration, whereas one cannot easily describe in geo- 
metric terms the scenario of an event horizon appearing and disappearing.) There is also no 
a priori assumption of equilibrium conditions: the concept of temperature is neither viable 
nor necessary, as is expected in all non-equilibrium conditions. Thermal or near-thermal ra- 
diance is a result of some specific conditions acting on the vacuum fluctuations in the system. 

4) We restrict our attention in this paper to near-thermal conditions because of tech- 
nical rather than conceptual limitations. In the near-thermal cases treated here, we want 
to add that the stochastic theoretical method is not the only way to derive these results. 
One can alternatively approach with the global geometric or thermal field methods, say, 
by working with generalized definitions of event horizons or quasi-periodic Green's func- 
tions. However, we find it logically more convincing and technically more rigorous to use the 
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stochastic method to define and analyze statistical concepts like fluctuations and dissipation, 
correlation and coherence. Certainly in the fully dynamical and non-equilibrium conditions, 
such as will be encountered in the full backreaction problem (not just confined to the linear 
response regime) this method is, in our opinion, more advantageous than the existing ones. 
Even though the technical problems will likely be grave, (because of the built-in balance 
between dissipation and fluctuations, as demanded by a self-consistent treatment), there are 
no intrinsic conceptual pitfalls or shortcomings. These issues and problems are currently 
under investigation. 
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